Discovering Strongly Correlated Quantum Spin Liquid 
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Strongly correlated Fermi systems are among the most intriguing and fundamental systems in 
physics. We show that the herbertsmithite ZnCu3(OH)6Cl2 can be viewed as a new type of strongly 
correlated electrical insulator that possesses properties of heavy-fermion metals with one exception: 
it resists the flow of electric charge. We demonstrate that herbertsmithite's low temperature proper- 
ties are defined by a strongly correlated quantum spin liquid made with such hypothetic particles as 
fermionic spinons which carry spin 1/2 and no charge. Our calculations of its thermodynamic and 
relaxation properties are in good agreement with recent experimental facts and allow us to reveal 
their scaling behavior which strongly resembles that observed in heavy-fermion metals. Analysis of 
the dynamic magnetic susceptibility of strongly correlated Fermi systems suggests that there exist 
at least two types of its scaling. 



PACS numbers: 75.40.Gb, 64.70. Tg, 76.60.Es, 71.10.Hf 

Strongly correlated Fermi systems represented by 
heavy-fermion (HF) metals arc well experimentally stud- 
ied but only recently have got an adequate theoretical 
description [l|. Landau Fermi liquid (LFL) theory is 
highly successful in the condensed matter physics. The 
key point of this theory is the existence of fermionic 
quasiparticlcs defining the thermodynamic, relaxation 
and dynamic properties of conventional metals. How- 
ever, strongly correlated Fermi systems encompassing a 
variety of systems that display behavior not easily un- 
derstood within the Fermi liquid theory and called non- 
Fermi liquid (NFL) behavior. A paradigmatic example of 
the NFL behavior is demonstrated by HF metals, where 
a quantum phase transition (QPT) induces a transition 
between LFL and NFL 0, 0]. QPT can be tuned by 
different parameters, such as the chemical composition, 
the pressure, and the magnetic field. Magnetic materials, 
in particular insulators, are interesting subjects of study 
due to a quantum spin liquid (QSL) that may develop in 
them, defining their low-temperature properties. Exotic 
QSL is formed with such hypothetic particles as fermionic 
spinons carrying spin 1/2 and no charge. A search for 
the materials is a challenge for condensed matter physics 
[|. In zero and high magnetic fields B [M3, the ex- 
perimental studies of herbertsmithite ZnCu3(OH)eCl2 
have discovered gaplcss excitations, analogous to quasi- 
particle excitations near the Fermi surface in HF met- 
als, indicating that ZnCu3(OH)gCl2 is the pro mising sys- 
tem to investigate its QPTs and QSLs |14Wl6| . The 
observed behavior of the thermodynamic properties of 
ZnCu3(OH)gCl2 strongly resembles that in HF metals 
since a simple kagome lattice has a dispersionless topolog- 
ically protected branch of the spectrum with zero excita- 
tion energy [l4, 17, 18 [. This indicates that QSL formed 
by the ideal kagome lattice and located near the fermion 



condensation quantum phase transition (FCQPT) can 
be considered as a strongly correlated quantum spin liq- 
uid (SCQSL). This observation allows us to establish a 
close connection between ZnCu3(OH)6Cl2 with its SC- 
QSL and HF metals whose HF system is located near 
FCQPT and, therefore, exhibiting an universal scaling 
behavior [J Q, 0. Thus, FCQPT represents QPT of 
ZnCu3(OH)gCl2 and both herbertsmithite and HF met- 
als can be treated in the same framework, while SCQSL 
is composed of spinons and these with zero charge and 
spin a = ±1/2 occupy the corresponding Fermi sphere 
with the Fermi momentum p F 

In our paper we show that both non-Fermi liquid and 
scaling behavior of such strongly correlated Fermi sys- 
tems as HF metals and ZnCu3(OH)eCl2 can be described 
within the frame of the theory of FCQPT. Analyzing ex- 
perimental data obtained in measurements on strongly 
correlated Fermi systems with different microscopic prop- 
erties we have found out that they demonstrate the uni- 
versal non-Fcrmi liquid behavior. Our analysis of the 
dynamic magnetic susceptibility of strongly correlated 
Fermi systems suggests that there exist at least two types 
of its scaling. We calculate the thermodynamic and relax- 
ation properties of herbertsmithite and HF metals. The 
calculations are in a good agreement with experimental 
data and allow us to detect the low-tempcraturc behavior 
of ZnCu 3 (OH) 6 Cl 2 defined by SCQSL as that observed 
in heavy fermion metals. 

To study theoretically the low temperature thermo- 
dynamic, relaxation and scaling properties of herbert- 
smithite, we use the model of homogeneous HF liquid 
[lj. This model permits to avoid complications associ- 
ated with the crystalline anisotropy of solids. Similar to 
the electronic liquid of HF metals, SCQSL is composed 
of chargeless fermions (spinons) with S = 1/2 occupy- 
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ing the corresponding two Fermi spheres with the Fermi 
momentum pp. The ground state energy E(n) is given 
by the Landau functional depending on the quasiparticle 
distribution function n CT (p), where p is the momentum 
and a is the spin index. The effective mass M* is gov- 
erned by the Landau equation P, Hij 
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Here we rewrite the quasiparticle distribution function 
as n*(p,B,T) = n a {p,B = 0,T = 0) + 5n*(p,B,T). 
The Landau amplitude F is completely defined by the 
fact that the system has to be at FCQPT [J [2(J-|22| , see 
20l - 22 1 for details of solving Eq. ([TJ . The sole role of Lan- 
dau amplitude is to bring the system to FCQPT point, 
where Fermi surface alters its topology so that the effec- 
tive mass acquires temperature and field dependencies. 
At this point, the term 1/M* vanishes, Eq. ([TJ becomes 
homogeneous and can be solved analytically [H [2(| . At 
B = 0, the effective mass, being strongly T - dependent, 
demonstrates the NFL behavior given by Eq. (QJ 



M*{T) ~ a T T 



-2/3 



(2) 



where ax is a constant. At finite T, under the applica- 
tion of magnetic field B the two Fermi spheres due to 
the Zeeman splitting are displaced by opposite amounts, 
the final chemical potential fi remaining the same within 
corrections of order B 2 . As a result, field B drives the 
system to LFL region, and again it follows from Eq. ([TJ 
that 



M*{B) ~ a B B 



-2/3 



(3) 



where as is a constant. It is seen from Eqs. ([2]) and ([3]) 
that effective mass diverges at FCQPT. At finite B and 
T, the solutions of Eq. (TJ M*(B,T) can be well ap- 
proximated by a simple universal interpolating function. 
This interpolation occurs between the LFL regime,given 
by Eq. © and NFL regime given by Eq. © (Jill. 
Experimental facts and calculations show that M*(B, T) 
as a function of T at fixed B reaches its maximum value 
M M at Tm, see e.g. [lj]. To study the universal scaling 
behavior of strongly correlated Fermi system, it is con- 
venient to introduce the normalized effective mass 
and the normalized temperature TV dividing the effec- 
tive mass M* and temperature T by their values at the 
maximum, M M and Tm respectively. In the same way, 
we can normalize other thermodynamic functions such as 
the spin susceptibility \ and the heat capacity C . As a 
result, we obtain 



X n (C/T) N ~ M* N , 



(4) 



where xn and (C/T)n are the normalized values of x 
and C/T, respectively. We note that our calculations of 



Mfl based on Eq. (QJ do not contain any fitting parame- 
ters. The normalized effective mass Mjy = M* /M M as a 

M 



function of the normalized temperature y = TV 
is given by the interpolating function [l[ 
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Here cq = (1 + C2)/(l + ci), c\ and C2 are fitting pa- 
rameters. Since magnetic field B enters Eq. OJ only in 
combination B[Ib /ksT ', we have T max oc B [l|,l20(, where 
fiB is the Bohr magneton and ks is the Boltzmann con- 
stant. Thus, for finite magnetic fields variable y becomes 



y = T/T/v oc k B T/n B B. 



(6) 



Since the variables T and B enter symmetrically Eq. ([5]) 
is valid for y = ^iBB/ksT. 

To construct the dynamic spin susceptibility 
X (q,^,T) = x'(q,w,T) + i X "(q,w,T) as a func- 
tion of momentum q, frequency uj and temperature T, 
again we use the model of homogeneous HF liquid lo- 
cated near FCQPT. To deal with the dynamic properties 
of Fermi systems, one can use the transport equation 
describing a slowly varying disturbance (5n CT (q, u>) of 
the quasiparticle distribution function no(p), and 
?i = 6 n + no- Wc consider the case when the disturbance 
is induced by the application of external magnetic field 
B = £> + \Bi(q,u) with B being a static field and 
XBi a w-dependent field with A — > 0. As long as the 
transferred energy uj < qpp/M* << u, where M* is 
the effective mass and /i is the chemical potential, the 
quasiparticle distribution function n(q, uj) satisfies the 
transport equation [23| 

dn 
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= qv p ^cr/j B (B + XBi) 



Here [Ib is the Bohr magneton and e v is the single- 
particle spectrum. In the field Bo, the two Fermi sur- 
faces are displaced by opposite amounts, ±Bo/Js, and 
the magnetization M. = ub(<5?i + — Sri—), where the two 
spin orientations with respect to the magnetic field are 
denoted by ±, and 5n± = ^2 p 5n±{p). The spin sus- 
ceptibility x is given by x = 9M/dB^ B=Ba . In fact, the 
transport equation ([7]) is reduced to two equations which 
can be solved for each direction ± and allows one to cal- 
culate 8n± and the magnetization. The response to the 
application of ASi(q, u>) can be found by expanding the 
solution of Eq. in a power series with respect to 
M*uj I (JPf- As a result, we obtain the imaginary part of 
the spin susceptibility 



x"(q,w) =/x| 



w(M*) 2 
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where F£ is the dimensionlcss spin antisymmetric quasi- 
particle interaction 1231. T he interaction Fq is found to 
saturate at Fg ~ -0.8 so that the factor (1 + F a ) 

in Eq. ((5]) is finite and positive. It is seen from Eq. © 
that the second term is an odd function of uj. There- 
fore, it does not contribute to the real part X ' an d forms 
the imaginary part x" ■ Taking into account that at rel- 
atively high frequencies u > qpp/M* <C fi in the hydro- 
dynamic approximation x' tx 1 /to 2 (2(| , we conclude that 
the equation 



x(q,w) 



M*p F 



7T 2 (1+F a ) 1+lTT 



q P F(l + Fg) 



(9) 



produces the simple approximation for the susceptibility 
X and satisfies the Kramcrs-Kronig relation connecting 
the real and imaginary parts of x- 

To understand how can x" an d X given by Eqs. ([5]) and 
(p?|), respectively, depend on temperature T and magnetic 
field B, we recall that near FCQPT point the effective 
mass M* depends on T and B. To elucidate a scaling 
behavior of x, we employ Eq. ([2]) to describe the tem- 
perature dependence of x- H follows from Eqs. (|9]) and 
^ that 



T 2/3 x(7» 



a i 



1 + ia 2 E 



(10) 



Here a\ and a 2 are constants absorbing irrelevant values 
and E = w/(/csT) 2 / 3 . As a result, the imaginary part 
x"(T,uj) satisfies the equation 



T 2 / 3 x"(7» 



a 3 E 



1 + a 4 E 2 ' 



[11) 



where (23 and a 4 are constants. It is seen from Eq. (| 1 1|) 
that T 2 / 3 x"(T,w) has a maximum (T 2 / 3 x"(T, w)) ma x 
at some E max and depends on the only variable -E. 
Equation ([TT1) is in accordance with the scaling be- 
havior of x"T 66 experimentally established in Ref. 
Q. As it was done for the effective mass when con- 
structing (|S|), we introduce the dimensionlcss function 
(T 2 ^ X ")n = T 2 / 3 x"/(T 2 / 3 X ")max and the dimension- 
less variable En ~ E/E max , and Eq. (fTTj) is modified to 
read 



(T 2 /V) 



b x E 



N 



N 



b 2 E 2 



(12) 



with b\ and b 2 are fitting parameters which are to adjust 
the function on the right hand side of Eq. (|12[) to reach 
its maximum value 1 at E n = 1. 

Now we construct the schematic T — B phase diagram 
of ZnCu 3 (OH) 6 Cl 2 reported in Fig. [TJ At T = and B = 
the system is near FCQPT without tuning. It can also 
be shifted from FCQPT by the application of magnetic 
field B. Magnetic field B and temperature T play the 
role of the control parameters, driving it from the NFL 
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FIG. 1: (color online). Schematic T — B phase diagram of 
ZnCu 3 (OH) 6 Cl2 located on the disordered side of FCQPT. 
The solid circle at the origin shown by the arrow represents 
FCQPT. Vertical and horizontal arrows show LFL-NFL and 
NFL-LFL transitions at fixed B and T respectively. The total 
width W of the NFL and the transition region is shown by the 
arrow. The inset demonstrates the behavior of the normalized 
effective mass versus normalized temperature Tjv as given 
by Eq. ([5J. Temperatures Tjv ~ 1 signify a transition region 
between the LFL regime with almost constant effective mass 
and NFL one, given by T~ 2 / 3 dependence. The transition 
region, where reaches its maximum at T/T max = 1, is 
shown by the arrows and hatched area both in the main panel 
and in the inset. 



to LFL regions as shown by the vertical and horizontal 
arrows. At fixed B and increasing T the system transits 
along the vertical arrow from the LFL region to NFL one 
crossing the transition region. On the contrary, at fixed T 
increasing B drives the system along the horizontal arrow 
from the NFL region to LFL one. The inset demonstrates 
the universal behavior of the normalized effective mass 
M^- versus normalized temperature Tjv as given by Eq. 
([5]). It follows from Eq. ([S]) and seen from Fig. [T]that 
the total width W of the NFL and the transition region, 
shown by the arrow in Fig. [TJ tends to zero at diminishing 
T and B since W oc T oc B. 

A few remarks are in order here. Equation (jlip is valid 
provided that the system approaches FCQPT from the 
disordered side as shown in the phase diagram [T] If the 
system is located on the ordered side then at B = the 
behavior of the effective mass as a function of T is given 
by 



M*(T) ~ a T T' 



(13) 



where a T is a constant. Upon taking into account Eq. 
(|13|) and acting in the same way as it was done when 
deriving Eq. (jllj). we obtain that the imaginary part 
x"{T,lv) is given by the equation 



T X "(T,oj) 



a 5 E 



E 2' 



(14) 
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FIG. 2: (color online). The experimental data on measure- 
ments of xn — (C/T)n — and our calculations of M% 
at fixed magnetic field are shown in the legends by points 
of different shape and solid curve respectively. It is clearly 
seen that the data collected on both ZnCu3(OH)6Ci2 [3 an d 
YbRli2Si2 merge into the same curve, obeying the scaling 
behavior. This demonstrates that SCQSL of herbertsmithite 
is close to FCQPT and behaves like HF liquid of YbRh 2 Si 2 
in magnetic fields. 



where CI5 and are constants, and E = ujjkBT. It is 
seen from Eq. (|14p that Tx"(T,uj) depends on the only 
variable E = uj/k B T. Thus, Eqs. dTTJ) and (14]) establish 
two types of scaling behavior of x"( w ) T). Since the scal- 
ing behavior of x"(w,T) is defined by the dependence of 
M* on T, one may expect new types of scaling especially 
at the transition region shown in Fig. [TJ 

Figure[3]reports the behavior of the normalized \n and 
specific heat (C/T)jv extracted from measurements on 
ZnCu 3 (OH) 6 Cl 2 and YbRh 2 Si 2 (27|, correspondingly. 
It follows from Fig. [2] that in accordance with Eq. ^ 
the behavior of xn coincides with that of (C/T)n in 
YbRh 2 Si 2 . 

In Fig. [3] the normalized xjv and (C/TU- extracted 
from measurements on ZnCu3(OH) 6 Cl 2 [3, Wi, are de- 
picted. To extract the specific heat C coming from the 
contribution of SCQSL from the total specific heat Ct{T) 
measured on ZnCu3(OH)fiCl 2 , we approximate Ct{T) at 



T > 2 K by the function 
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C t (T) = ai T 3 + a 2 T 



1/3 



(15) 



where the first term proportional to a\ is due to the lat- 
tice (phonon) contribution and the second one is deter- 
mined by SCQSL when it exhibits the NFL behavior as 
it follows from Eq. C oc TM* cx T 1 / 3 . Taking into 
account that the phonon contribution is B-field indepen- 
dent, we obtain C(B,T) = C t (B,T) - aiT 3 . It is seen 
from Figs. f2]and[3]that (C/T)jv — Xn displays the same 



FIG. 3: (color online). The normalized susceptibility xn — 
M^f and the normalized specific heat (C/T)n — M"^ of SC- 
QSL versus normalized temperature TV as a function of the 
magnetic fields shown in the legends, xn and (C/T)n are 
extracted from the data of [3] and [H, respectively. 



scaling behavior as {C/T)m measured on the HF metal 
YbRh 2 Si 2 . Therefore, the scaling behavior of the ther- 
modynamic functions of herbertsmithite is the intrinsic 
feature of the compound and has nothing to do with mag- 
netic impurities. The observed scaling behavior of xn 
and (C/T)jv in magnetic fields shown in Figs. [2] and [3] 
rules out a possible supposition that extra Cu spins out- 
side the kagome planes considered as weakly interacting 
impurities could be responsible for the divergent behavior 
of the susceptibility at low temperatures. 

In Fig. |4l consistent with Eq. (fL2|) . the scaling 
of the normalized dynamic susceptibility (T 2 / 3 x")n ex- 
tracted from the inelastic neutron scattering spectrum 
of both herbertsmithite 0] , Panel A, and the HF metal 
Ceo.g 2 5Lao.o75Ru 2 Si 2 [2§|, Panel B, is displayed. In Fig. 
31 Panel C, the dynamic susceptibility (Tx") extracted 
from measurements of the inelastic neutron scattering 
spectrum on the HF metal YbRh 2 Si 2 [3(j is shown. 
The data {Tx") exhibit the scaling behavior over three 
decades of the variation of both the function and the 
variable, thus confirming the validity of Eq. (fT4]) . The 
scaled data obtained in measurements on such quite dif- 
ferent strongly correlated systems as ZnCu3(OH)eCl 2 , 
Ceo.g 2 5Lao.o75Ru 2 Si 2 and YbRh 2 Si 2 collapse fairly well 
onto a single curve over almost three decades of the scaled 
variables. It is seen that our calculations shown by the 
solid curves are in good agreement with the experimental 
facts. Some remarks on a role of both the disorder and 
the anisotropy are in order. The anisotropy is supposed 
to be related to the Dzyaloshinskii-Moriya interaction, 
exchange anisotropy, or out-of-planc impurities. Mea- 



5 



ZnCu (OH) CI 



■ 


T=77mK 


o 


T=280mK 


A 


T=1K 


▼ 


T=3.5K 


< 


T=12K 


► 


T=42K 




1.0 



0.5 



Q, =(0.6910) 



o 


T=5K 


A 


T=12.5K 


O 


T=16,5K 


□ 


T=24K 




T=35K 





T=53K 


□ 


T=80K 





FIG. 4: (color online). Panels A and B, (T 2/3 x")n plot- 
ted against the unitless ratio En = w/((fcsr) 2//3 i5 max ). The 
data are extracted from measurements on ZnCu3(OH)eCl2 
(tJ, Panel A, and Panel B, on Ceo.925Lao.o75Ru 2 Si2 obtained 
at Qi HI]. Panel C, T X " plotted against E = oj/k B T. The 
data are extracted from measurements on YbRh2Si2 [30|. The 
solid curves, Panels A and B, are fits with the function given 
by Eq. (O, Panels C, with the function given by Eq. JT3}. 



surements of the susceptibility on the single crystal of 



herbertsmithite have shown that it closely follows that 
measured on a powder sample [H, At low temper- 
atures T < 70 K, the single-crystal data do not show 
substantial magnetic anisotropy 0, These confirm 
that the stoichiometry, disorder and anisotropy do not 
contribute significantly to the results at relatively low 
temperatures. As we have seen above, the scaling be- 
havior of the thermodynamic and relaxation properties 
of herbertsmithite is its intrinsic feature and has noth- 
ing to do with the impurities 15] . These observations 



are in agreement with ageneral consideration of scaling 
behavior of HF metals [l[ . 

In summary, we have considered the non-Fermi liquid 
behavior and the scaling one of such strongly correlated 
Fermi systems as insulator ZnCu3(OH)6Cl2 and HF met- 
als Ceo.925Lao.o75R-U2Si2 and YbRti2Si2, and shown that 
these are described within the frame of the theory of 
FCQPT. Our analysis of the dynamic magnetic suscepti- 
bility of strongly correlated Fermi systems suggests that 
there exist at least two types of its scaling. We calculate 
the thermodynamic and relaxation properties of herbert- 
smithite and HF metals. The calculations arc in a good 
agreement with experimental data and allow us to iden- 
tify the low-temperatures behavior of ZnCu3(OH) 6 Cl2 
determined by SCQSL as that observed in heavy fermion 
metals. Thus, herbertsmithite can be viewed as a new 
type of strongly correlated electrical insulator that pos- 
sesses properties of heavy-fermion metals with one ex- 
ception: it resists the flow of electric charge. 
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